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I .  INTRODUCTION 

fame  years  ago,  a  number  of  people  (1-4)  Independently 
suggested  the  use  of  a  saturable  absorber  to  induce  optical  bistabili¬ 
ty  (hereafter  referred  to  as  OB).  The  phenomenon  was  demonstrated  and 
studied  recently  by  Gibbs,  McCall  and  Venkatesan  (5)  who  used  a  cw  dye 
laser  to  excite  atoms  of  sodium  vapor  in  a  cell  between  the  plates  of 
a  Fabry  Perot  interferometer.  Under  certain  experimental  and  material 
conditions,  they  observed  a  nonlinear  dependency  of  the  transmitted 
field  as  a  function  of  the  dye-laser  input  which  exhibited  hysteresis, 
differential  gain,  and  bistable  behavior. 

Interest  in  OB  from  the  practical  applications  standpoint 
stems  from  obvious  device  applications  as  the  optical  analog  of  the 
transistor,  optical  clipper  or  limiter  and  digital  optical  memory  ele¬ 
ment.  The  recent  demonstration  of  the  phenomenon  (5)  has  generated 
considerable  interest  from  the  fundamental  point  of  view  as  a  clear 
example  of  spontaneous  ordering  in  an  open,  stationary  system  of  mat¬ 
ter  interacting  with  light.  This  has  led  to  a  surge  of  recent  theo¬ 
retical  activity\^/''\^ 

The  first  model  for  OB  was  proposed  by  McCall  (6) ,  who  in¬ 
troduced  a  nonlinear  susceptibility  into  the  Maxwell-Bloch  represen¬ 
tation  in  the  full  propagation  treatment.  His  results  are  necessarily 
calculatlonal  rather  than  analytical.  Bonifacio  and  Lugiato  (7)  were 
the  first  to  obtain  the  main  features  of  OB  from  Implementation  of  the 
Maxwell-Bloch  model,  which  they  solved  analytically  in  closed  form. 
They  later  extended  their  model  (8)  to  a  quantum-statistical  represen¬ 
tation  of  a  stationary  system  far  removed  from  thermodynamic  equilib¬ 
rium,  but  with  the  same  essential  assumptions,  the  most  crucial  being 
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the  "mean  field"  approximation  which  amounts  to  requiring  that  the 
field  be  sufficiently  uniform  over  the  volume  of  active  atoms.  Others, 
such  as  Walls  et  al.  (9)  have  considered  similar  models  which  are,  on 
the  whole,  in  agreement  with  the  results  of  Bonifacio  and  Lugiato. 

All  of  these  models  have  one  essential  comnon  characteris¬ 
tic,  namely,  they  consider  OB  within  the  framework  of  individual  or 
Independent  atom  response  to  the  electromagnetic  field.  These  are 
therefore  necessarily  absorptive  models  which  describe  OB  as  a  tran¬ 
sition  from  a  state  of  coherent  response  of  individual  atoms  to  the 
externally- applied  driving  field  and  high  absorption  (low  transmis¬ 
sion)  to  a  state  of  saturation  or  bleaching  of  the  material  (high 
transmission) .  The  experimental  results  on  sodium  vapor  (5)  exhib¬ 
ited  a  large  dispersive  contribution,  and  these  models  do  not  appear 
to  describe  the  dispersive  aspects  even  qualitatively,  although  the 
absorptive  characteristics  appear  to  be  adequately  qualitatively  de¬ 
scribed. 


It  seemed  to  us  that  OB  should  stem,  at  least  under  certain 
conditions,  from  the  same  kind  of  fundamental  interaction  that  gives 
rise  to  the  dynamical  cooperative  process  known  as  superfluorescence 
(10,  11).  This  prompted  us  to  take  a  fresh  approach  to  the  interpre¬ 
tation  of  OB  in  terms  of  a  model  which  emphasizes  the  contribution 
from  interatomic  cooperative  processes  in  matter  interacting  with 
light.  Our  model  (12)  differs  from  other  previous  models  in  that  it 
Includes  the  possibility  for  atom-atom  interaction  and  correlation 
via  the  mutual  Internal  electromagnetic  field  (virtual  photon  ex¬ 
change)  as  well  as  the  individual  atom  coupling  to  the  electromagnet¬ 
ic  field. 

The  results  of  our  model  (12)  and  other  models  (7-9)  pre¬ 
dict  identical  behavior  for  "absorptive"  OB  (perfect  cavity,  laser 
and  atomic  tuning).  However,  our  results  fit  the  dispersive  char¬ 
acteristics  reported  for  the  sodium  experiment  (5)  qualitatively 
quite  well.  In  addition,  we  predict  new  results  in  reversible  hys¬ 
teresis  and  bistability  in  the  variation  of  Inverse  effective  temper¬ 
ature  Ts  with  the  output  field  for  fixed  input  above  a  certain  thres¬ 
hold  value.  This  new  result  (13)  has  a  wide  range  of  potential 
applications,  including  a  sensitive  IR,  FIR  or  millimeter  wavelength 
detector  with  frequency  upconverslon. 

In  addition,  our  results  provide  an  interpretation  of  OB  as 
a  true  first-order  phase  transition,  from  atomic  order  to  disorder 
due  to  virtual  photon  exchange  between  atoms,  far  from  thermodynamic 
equilibrium.  In  the  absence  of  an  applied  field,  our  results  show 
the  existence  of  a  second-order  phase  transition  (12,14)  below  a 
critical  temperature  determined  by  material  parameters  and  cavity 
detuning. 
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Another  unique  characteristic  of  our  results  is  that  the 
model  constitutes  the  optical  analog  of  the  Meisner  effect  (15)  in 
superconductivity.  Due  to  the  collective  atomic  behavior  (optical 
Meisner  effect) ,  our  results  predict  for  the  first  time  the  existence 
of  OB  in  a  small  volume  without  mirrors.  This  could  have  far-reach¬ 
ing  applications  in  the  use  of  arrays  of  micron  size  elements  for 
optical  digital  information  storage  and  retrieval  as  well  as  optical 
digital  image  formation. 

Finally,  it  should  be  pointed  out  that  this  work  is  a  di¬ 
rect  outgrowth  of  our  program  in  swept-gain  superradiance  (16-18). 
Furthermore,  Ml COM  will  co-host,  with  ARO-D,  an  International  Con¬ 
ference  and  Workshop  on  Optical  Bistability,  in  Asheville,  North 
Carolina,  3-5  June  1980.  This  conference  is  expected  to  significant¬ 
ly  advance  the  field  and  to  provide  a  timely  record  of  the  state-of- 
the-art  in  this  rapidly-developing  and  potentially  extremely  useful 
field.  The  conference  and  its  impact  are  anticipated  to  be  similar 
and  of  much  the  same  significance  as  the  conference  hosted  by  MICOM 
on  Superradiance  (19) . 

II.  THEORETICAL  MODEL  AND  INTERPRETATION  OF  PHENOMENON 

We  consider  a  collection  of  N  identical  two-level  atoms 
with  energy  level  separation  e  interacting  with  the  electromagnet¬ 
ic  field  in  the  dipole  approximation.  The  electromagnetic  field 
which  Interacts  with  the  atoms  is  comprised  of  an  internal  field  and, 
in  addition,  an  applied  cw  field  assumed  to  be  in  a  coherent  state 
(20),  with  field  amplitude  a  and  carrier  frequency  uo>  The  atoms 
of  the  material  are  considered  to  be  confined  to  a  volume  Vg  which  is 
much  smaller  than  the  resonance  wavelength  Xr  associated  with  the 
atomic  bare-state  separation.  This  assumption  does  not  critically 
affect  our  results  and  tends  to  allow  emphasis  upon  the  nature  of 
collective  atomic  aspects  in  the  results  in  an  unencumbered  fashion. 

The  Hamiltonian  H  which  represents  this  model  is  written 
in  the  form 

H  -  H 

o 

where,  in  units  such  that  ft  ■  c 

ho  ■  £  vl\  * 

k*o 


+  H’  (II-l) 

-  1  , 

N 

*se  £  aj  (II-2a) 

4ml  J 
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(N®  \  /  lilt  N  \ 

£  E  ^<5  +  k-cj+|a*e  °  £  a"  +  h  •  c  ,  (II-2b) 

where  h  •  c  means  hermetian  conjugate,  and  we  have  used  the  rotating 
wave  approximation. 

The  first  and  second  terms  in  (II-2a)  represent  the  free 
field  and  atomic  system,  respectively,  whereas  the  first  two  terms 
of  (II-2b)  describe  the  interaction  of  the  atoms  with  the  internal 
field  and  the  last  two  terms  describe  their  interaction  with  the  ex¬ 
ternally-applied  coherent  field.  The  factor  a  is  the  complex  extern¬ 
ally-applied  field  amplitude.  The  field  photon  creation  and  annihil¬ 
ation  operators  conform  to  the  Bose  conmutation  relations  (12),  and 
the  atomic  operators  are  the  usual  Pauli  operators,  and 

gk  -  -i  (Wj^V)1*  ^  •  d  (H-3) 

where  V  is  the  quantization  volume  and  and  d^  are  the  polarization 
vector  and  the  matrix  element  for  the  transition  dipole  moment  for 
the  atoms  (21) . 

We  are  Interested  in  determining  the  steady  state  properties 
of  the  system  described  by  (II-2a,  ll-2b).  In  order  to  retain  the 
essential  many-body  aspects  of  the  system,  it  is  convenient  to  treat 
the  system  in  an  ensemble  representation.  This  means  taking  ensemble 
averages  of  relevant  dynamical  operators  such  as  population  differ¬ 
ence  and  dipole  moment  operators.  To  facilitate  this,  we  perform  the 
following  operations:  (i)  Canonically  transform  to  a  perpendicularly 
doubly  rotating  frame,  rotating  at  the  laser  frequency  ,  so  the 
Hamiltonian  (II-2)  becomes  explicitly  time- independent  (12) .  The 
electromagnetic  field  then  appears  as  a  d.c.  field  in  the  transformed 
frame,  (ii)  Introduce  an  effective  temperature  Ts  in  the  transform¬ 
ed  representation  which  is  the  exact  analog  of  a  "spin  temperature" 
(12,22).  This  is  equivalent  to  a  definite  statement  of  the  form  of 
the  density  operator  p  ,  in  the  transformed  representation. 

Although  the  system  is  an  open  system  far  removed  from 
thermodynamic  equilibrium,  we  can  analyze  the  macroscopic  properties 
of  the  system  in  the  transformed  representation  in  terms  of  the  free 
energy  F,  given  by  _g  F  H 

e  8  -  tr  e  8  T  ,  (II-4) 

where  8g  ■  l/kTg  and  k  is  Boltzmann's  constant. 
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Making  use  of  the  transformation  to  the  doubly  rotating 
frame  (2-4)  of  reference  (12)  and  adiabatically  eliminating  the  in¬ 
ternal  field  variables,  the  transformed  Hamiltonian  H-p  ,  through 
second  order  in  the  atom-field  interaction,  was  shown  by  Bowden  and 


Sung  (12) 

to 

be 

*T-*oT  + 

H' 

T 

(II-5) 

where 

N 

HqT  «  S2a+a  + 

^a  £ 
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■J  ■ 

(II-6a) 

N  4. 

N 
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Q 
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(II-6b) 

j  *  l 

Here, 

Q 

» 

e 

oE 

«• 

fia  " 

e  -  a) 
o 

• 

(XI-7) 

For  simplicity,  we  have  considered  only  a  single  radiation  field  mode 
of  frequency  <u  Interacting  with  the  atomic  transition.  If  the  atomic 
medium  is  in  a  cavity,  then  u  is  the  cavity  mode  frequency.  If  there 
is  no  cavity,  then  u  -  £  .  The  first  term  in  (II-6b)  is  interpreted 
as  atomic  interaction  via  the  electromagnetic  field,  where  the  effec¬ 
tive  atom-atom  coupling,  g  ,  is  field-dependent  (12), 


where 


In  the  above. 


8  - 


(<|3— (i)  -  X) 

o 

(to-w — X  )^  +  Y^ 
o 


(II-8) 


go  "  ^Td"k  P  *  <H-9) 

X  is  the  Stark  shifted  atomic  level  separation  (12), 


(11-10) 


where  |&|  is  the  Stark  shift  which  is  dependent  upon  the  value  for 
the  Internal  field  intensity.  Since  the  latter  is  to  be  determined 
in  the  result,  ] A ]  must  be  determined  self-conslstently  later  on  in 
the  results  of  the  calculation.  It  is  crucial  to  retain  the  Stark 
shift  in  (11-10)  (usually  neglected  in  second-order  calculations)  in 
order  to  obtain  the  first-order  phase  transition  which  leads  to  OB. 
The  damping  term  Y  appearing  In  (II-8)  is  just  the  cavity  width 
(12)  or,  if  no  cavity  is  present,  it  is  just  the  photon  escape  rate 
from  the  atomic  volume.  In  (II-9),  d  is  the  matrix  element  of  the 
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transition  dipole  moment,  k  is  the  density  of  field  modes  in  the 
neighborhood  of  the  principal  mode  of  frequency  u  ,  and  p  is  the 
total  number  of  atoms  N  in  a  cubic  wavelength  at  atomic  resonance  Xr. 

The  equation  (II-6)  is  our  working  Hamiltonian.  The  atom- 
radiation  field  Interaction  (ll-2b)  has  now  been  recast  in  terms  of  an 
effective  interatomic  pair  interaction  via  the  internal  field  (II-6b). 
The  pair  correlation  is  depicted  by  the  first  term  in  (II-6b)  and  may 
lead  to  cooperative  atomic  effects  by  virtual  photon  exchange.  This 
term  is  identical  to  the  interaction  which  leads  to  cooperative  atom¬ 
ic  behavior  in  superfluorescence,  (10,23).  The  inclusion  of  the  Stark 
shift  correction  in  g  ,  (II-8) ,  (II- 10)  gives  negligible  contributions 
to  the  results  for  superfluorescence,  but  is  absolutely  essential  to 
arrive  at  OB  in  the  steady  state  for  this  model.  This  causes  the 
effective  interatomic  coupling  to  be  internal  field  intensity  depen¬ 
dent,  and  thus  leads  to  feedback  between  the  atoms  and  their  mutual 
radiation  field.  The  atomic  pair  correlation  term  in  the  interaction 
is  also  Identical  to  that  for  Cooper  pairs  in  the  BCS  theory  of  super¬ 
conductivity  (12,15). 

Ue  proceed  now  to  calculate  the  free  energy  F  from  (II-4). 
To  perform  the  trace  on  the  right-hand  side  of  (II-4),  it  is  neces¬ 
sary  to  linearize  the  Hamiltonian  (II-6)  and  to  obtain  a  thermodynam¬ 
ically  equivalent  linearized  Hamiltonian  (25,26).  The  method  of  ref¬ 
erences  (25)  and  (26)  is  used  to  obtain 

HJ.  -  N  I  jv|2  -  N  A*  -  NA<f  ,  (11-11) 

where 

A  =  i  V*  -  a*  .  (11-12) 

Here,  v  is  a  variational  parameter  used  in  the  linearization,  to  be 
determined  by  minimization  of  the  free  energy  F  ,  (25,26). 

In  the  above  expressions,  we  have  used  the  assumption  that 
all  the  atoms  are  Identical  to  drop  the  subscripts,  j  ,  and  to  re¬ 
place  the  sums  by  N  and  the  double  sums  by  N(N-l). 

The  thermodynamically  equivalent  Hamiltonian  (11-11) ,  l.e., 
the  linearized  Hamiltonian  which  yields  the  same  free  energy  F  in  the 
limit  of  large  N  as  the  original  Hamiltonian  (II-6b)  has  an  imsiediate 
physical  interpretation.  The  interaction  expressed  in  (II-6b)  is 
thermodynamically  equivalent  to  a  single  atom,  mean  field  interaction 
expressed  in  (11-11) .  The  factor  A  ,  (11-12),  is  the  mean  field,  l.e., 
the  total  field  acting  on  a  single  atom  due  to  the  combined  atomic  re¬ 
action  fields  of  all  the  other  atoms  (first  term  in  (11-12)  and  the 
local  value  of  the  externally-applied  field  (second  term  in  (11-12)). 
The  expression  for  the  total  local  field  at  an  atomic  site, (11-12) , 
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is  just  Maxwell's  equation  (all  electromagnetic  fields  appear  as  d. c. 
fields  in  the  transformed  frame  rotating  at  the  laser  frequency  u>0) . 

The  free  energy  F  is  now  easily  evaluated  from  (II— 4)  us¬ 
ing  (II- 11)  to  give. 


—  2  W  B_  i  ~ B 

F  -  N  g  |v|  -  log  [2  cosh  ~y-\1+  -j-  log  (1-e  8  ) 


(11-13) 


The  variational  parameter  v  is  evaluated  from  the  normal  equation, 
(25,26), 

3F/3v  -  0  . 


If  we  let 


x  “  A 


equation  (11-14)  yields 


y  •  -  a 
B 


3F  *  _  *  tan^  2  * 

£  “  "IV  -  2  g  *  - ^ -  -  0 


3v 

From  (11-16), 


v  -  2  x 


tanh  ~  X 


(11-14) 

(H-15) 

(11-16) 

(11-17) 


which  is  easily  shown  to  be  the  ensemble  average  fox  the  macroscopic 
atomic  polarization  (12).  It  is  furthermore  easily  shown  that  the 
above  condition  obtained  from  (11-14)  corresponds  to  an  absolute 
minimum,  l.e., 

32F/3vZ  >  0  .  (11-18) 

It  is  to  be  noted  that  the  Stark  shift  A  is  now  determined  self- 
consistently  from  (11-12)  and  (11-17). 

If  we  use  the  definitions  (11-15)  in  (11-12)  and  use 
(11-12)  to  eliminate  V  from  (11-17),  then  after  multiplying  both 
sides  of  (11-17)  by  g,  the  following  equation  of  state  for  the  appl¬ 
ied  field  y  as  a  function  of  the  internal  field  x  is  obtained, 

tanh  £*■  X 

x  -  y  -  2  g  x  - j -  .  (11-19) 


This  is  the  main  result  of  this  paper  and  leads  to  bistability,  re¬ 
versible  hysteresis  and  differential  gain  between  the  three  quanti¬ 
ties,  x,  y,  and  B  .  These  characteristics  will  be  discussed  in  the 
next  section. 
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III.  INTERPRETATION  OF  RESULTS 
A.  Optical  Bistability  as  a  First-order  Phase  Transition 

The  main  result  of  our  model,  (11-19),  has  the  form  typical  of  a 
first-order  phase  transition.  For  absolute  zero  of  effective  tem¬ 
perature  and  for  perfect  laser  and  cavity  tuning,  i.e.,  for  Tg  0  , 
to  *  e  ■  u)0  ,  (11-19)  reduces  to  exactly  the  form  for  the  equation 
of  state  obtained  earlier  from  the  single- atom,  coupled  Maxwell- 
Bloch  models  (8), (9),  i.e., 

y  -  x  +  2c - - - s-  ,  (III-l) 

(1  +  x)Z 

where 

c  -  g^/y2  ,  (HI-2) 

and  we  have  written  x  and  y  in  units  of  y  .  Here,  however,  the  phys¬ 
ics  is  quite  different  as  will  be  discussed  later. 

Equation  (11-19)  is  shown  plotted  in  Figure  1,  input  field 
amplitude  y  vs.  output  field  amplitude  x  for  various  values  of  equal 
cavity  and  laser  offtuning.  It  is  readily  established  that  the  mini¬ 
mal  condition  for  OB  to  occur  is  that  c  >  4,  (8,12). 

It  is  easily  shown  (27)  that  all  points  on  the  curves  in 
Figure  1  are  stable  steady-state  conditions  except  those  for  which 
dx/dy  <  0.  Thus,  the  points  on  the  curves  between  the  turning  points 
having  negative  slope  are  unstable.  Increasing  the  applied  field  y 
from  zero,  the  output  follows  a  nearly  linear  response  until  the 
first  turning  point  is  reached  along  the  curve  for  the  equation  of 
state.  In  this  regime,  the  material  responds  mainly  hy  atom-atom 
coupling  via  virtual  photon  exchange,  i.e.,  g  ,  (II-8)  is  large  due 
to  the  relatively  small  Stark  shift  (11-10) .  The  reaction  field  due 
to  the  collective  response  of  the  atoms  (first  term  on  the  right-hand 
side  of  (11-12))  opposes  the  applied  field  (second  term  on  the  right- 
hand  side  of  (11-12))  to  produce  a  small  output  field  (11-12,(11-15). 
This  phenomenon,  due  to  collective  response  to  an  external  field,  is 
somewhat  analogous  to  the  Meisner  effect  (28)  in  the  theory  of  super¬ 
conductivity.  In  this  regime  the  macroscopic  transverse  dipole  mom¬ 
ent,  V  ,  (11-17)  is  small,  again  due  to  the  relatively  small  Stark 
shift. 

Once  the  first  turning  point  is  reached  (point  A  of  curve 
a),  by  increasing  y,  there  are  no  steady  state  values  available  for 
the  system  if  y  is  further  increased,  unless  the  system  jumps  to  a 
higher  value  of  the  output  (point  B  of  curve  a).  This  jump 
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corresponds  to  a  sudden  decrease  in  g  ,  the  atom-atom  coupling,  and 
a  sudden  increase  in  the  macroscopic  polarization, v  .  Further  in¬ 
crease  in  y  produces  a  nearly  linear,  empty  cavity  response  In  x. 
That  is ,  the  material  is  nearly  transparent  to  the  input  y ,  i . e . , 
x  z  y.  In  this  regime,  the  atom-atom  pair  correlation  is  destroyed 
and  each  individual  atom's  dipole  moment  follows  the  driving  field 
y,  coherently  with  definite  phase.  This  could  be  called  the  regime 
of  "cw  self- induced  transparency". 


Fig.  1,  Effect  of  cavity  and 
laser  detuning  on  optical  bi¬ 
stability.  Input  field  y  vs. 
internal  field  x  from  (11-19). 

Curve  a:  c  *  40,  6m  V  m  0; 

curve  b:  c  ■  40,  <5*  v**  2; 

curve  c:  c  ■  40,  6  ■*  v  “  4. 


6  -  (w-u)0)/y;  v-  (e-oj0)/y; 
T  -  0. 


8 


Fig.  2.  Asymmetric  effect  of 
cavity  detuning  on  optical  bi¬ 
stability.  Input  field  y  vs.  in¬ 
ternal  field  x  from  (11-19) . 
Curve  a:  c  ■  60,  10;  curve 

b:  same  parameters  as  for  curve 
a,  except  6 ■  -9;  curve  c:  same 
parameters  as  for  curve  a,  ex¬ 
cept  6  *  9;  curve  d:  same  param¬ 
eters  as  for  curve  a,  except 

6  -  8;  T  -  0. 

8 


The  asymmetric  effects  of  cavity  offtuning  for  fixed  laser 
frequency  is  shown  in  Figure  2.  These  results  exhibit  large  dis¬ 
persive  contributions  for  cavity  offtuning  and  are  in  good  qualita¬ 
tive  agreement  with  the  experimental  results  for  sodium  vapor  (5). 
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The  effective  temperature  dependence  of  the  hysteresis  cy¬ 
cle  and  OB  threshold  is  shown  in  Figure  3  for  input  y  vs.  output  x 
from  (11-19)  in  units  of  Y  .  For  convenience,  we  have  taken  perfect 
tuning,  so  (11-19)  becomes,  in  units  of  Y  , 


y  *  x  +  2c - — — a—  tanh  \  Tx  ,  (III-3) 

(14*) 


where 

r  -  YBS  =  Y/kTg  .  (III-4) 

B.  Effective  Temperature-induced  Bistability  and  Hysteresis  and  the 
Optical  Meisner  Effect 


It  is  interesting  to  consider  the  inversion  of  (III-3)  to  obtain 
T  in  terms  of  the  input  y  and  output  x,  (29), 

r  -inaffr*.*  >  .  (III-5) 

x  (2cx  -  (y-x) (x  +1)  I 

Equation  (III-5)  is  plotted  in  Figure  4  for  the  inverse  temperature, 
r  vs.  output  field  x  for  fixed  y  in  units  of  y  .  Each  curve  corres¬ 
ponds  to  a  different  value  for  the  fixed  input  y  *  yc.  It  is  seen 
that  bistability  and  hysteresis  occur  in  the  output  field  x  by  vary¬ 
ing  the  inverse  temperature  F  when  the  input  y  is  fixed  at  a  value 
above  a  certain  threshold. 

Consider  the  curve  in  Figure  4  showing  the  largest  hyster¬ 
esis  area,  corresponding  to  the  appropriate  fixed  input  field  y  =  yc. 

As  the  effective  temperature  Ts  is  increased  (  T  decreased)  from  some 
small  value  to  larger  values,  the  output  x  increases  slightly  as  the 
reaction  field  gv  in  (11-12),  due  to  the  collective  atomic  pair  in¬ 
teraction  which  opposes  the  applied  field  a  ,  is  diminished  as  the 
temperature  Ts  is  increased,  (11-17).  This  is  due  to  simultaneous 
reduction  of  the  macroscopic  dipole  moment  (11-17)  because  of  the  in¬ 
crease  in  effective  temperature  Ts  ,  and  to  the  decrease  in  g  (II-8) 
as  a  consequence  of  the  increase  in  the  field-induced  Stark  shift, 

(11-10),  (11-15). 

As  T  is  decreased  to  the  value  corresponding  to  point  A 
in  Figure  4,  any  further  decrease  in  T  will  necessarily  correspond  to 
a  stationary  state  on  the  upper  part  of  the  curve,  point  B.  The  sys¬ 
tem  undergoes  an  abrupt  increase  in  the  output  x  as  a  consequence  of 
the  catastrophic  breakdown  of  the  atomic  pair  correlation  and,  hence, 
the  reaction  field  gv  (11-12).  In  other  words,  along  the  path  AB, 

Figure  4,  the  macroscopic  transverse  polarization  V  (11-17)  increases 
abruptly  due  to  the  increased  Stark  shift  caused  by  the  transition  to 

i 

i 
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larger  total  local  field.  In  the  same  process,  however,  the  atomic 
pair  correlation  g  undergoes  a  catastrophic  reduction  in  the  trans¬ 
ition  from  A  to  B.  As  T  is  further  decreased,  the  macroscopic  pol¬ 
arization  diminishes  due  to  the  increase  in  Ts  and  the  system 
approaches  the  state  of  atomic  saturation  (bleaching)  at  high  abso¬ 
lute  temperature.  The  reverse  process  takes  place  as  T  is  in¬ 
creased  (Ts  decreased). 


Fig.  3.  Normalized  externally- 
applied  field  y  vs.  normalized 
internal  field  x  according  to 
Eq.  (III-3)  with  c  *  60.  Each 
curve  corresponds  to  a  diffe¬ 
rent  value  for  T,  and  in  order 
of  decreasing  values  for  the 
threshold  value  for  y,  T-1  1.00, 
0.65,  0.35,  0.25.  The  threshold 
value  of  T  for  bistability  in 
this  case  is  T  *  0.35.  The  dot¬ 
ted  lines  and  arrows  indicate 
a  particular  hysteresis  cycle. 


Fig.  4.  Tvs.  x/yc  according  to 
Eq.  (III-5)  for  c  =  60  and  var¬ 
ious  values  for  the  fixed  exter¬ 
nally-applied  field  yc.  In  order 
of  decreasing  values  for  the 
threshold  for  T,  yc  “  21.8,  21.7, 
21.6,  21.3,  21.0,  20.0.  The  thres¬ 
hold  value  of  yc  for  bistability 
is  yc  =  21.0  which  is  seen  to  be 
consistent  with  the  threshold 
value  of  T  for  bistability  in 
Fig.  3.  The  dotted  lines  and 
arrows  indicate  a  particular 
hysteresis  cycle. 
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The  reduction  of  the  internal  field  x  due  to  atomic  pair 
correlation  (11-12) ,  and  the  abrupt  reduction  of  the  reaction  field 
at  a  critical  temperature,  Ts,  with  consequent  increase  in  the  total 
internal  field  x,  for  fixed  input  yc  shown  in  Figure  4  is  quite  anal¬ 
ogous  to  the  Meisner  effect  in  the  theory  of  superconductivity  (15,28). 
We  are  tempted  to  label  this  effect,  the  "optical  Meisner  effect". 

C.  Second-order  Phase  Transition  in  the  Absence  of  Applied  Field 

In  the  absence  of  applied  field,  i.e.,  for  y  =  0,  (11-19)  pre¬ 
dicts  the  existence  of  a  second-order  phase  transition,  provided  that 
the  following  gap  equation  is  satisfied  (12), 

tanh  *a0  e  =  -y-  »  (HI-6) 

2go 

where  Tc  is  here  the  actual  temperature  of  the  heat. bath. 

If  this  equation  can  be  satisfied,  it  defines  a  critical 
temperature  Tc  below  which  the  system  exhibits  a  macroscopic  polari¬ 
zation  (ordered  phase)  and  above  which  there  is  no  macroscopic  pol¬ 
arization  (disordered  phase) .  It  is  readily  seen  that  (III-6)  can¬ 
not  be  satisfied  without  off tuning  provided  by  cavity  mirrors. 
Furthermore,  it  is  readily  determined  that  for  conditions  where  OB  is 
observed  in  the  optical  regime,  i.e.,  c  >  4,  (III-2) ,  (III-6)  will 
never  be  satisfied  (12).  However,  it  can  be  satisfied  for  microwave 
frequencies.  This  explains  why  this  effect  has  not  yet  been  observed 
in  OB  experiments. 


IV.  SUMMARY  AND  CONCLUSIONS 

The  effort  at  MICOM  in  quantum  optics  has  produced  results 
which  for  the  first  time  predict  OB  which  is  caused  by  cooperative 
atom-atom  interaction  via  their  mutual  radiation  fields  (12) .  The 
mechanism  for  the  atomic  pair  correlation  can  be  interpreted  as  due 
to  virtual  photon  exchange  interaction  in  analogy  with  the  inter¬ 
action  between  Cooper  pairs  via  the  phonon  field  in  the  BCS  theory  of 
superconductivity  (30).  The  measure  of  the  strength  of  the  atomic 
pair  correlation  is  g  (II-8).  The  main  result  of  our  theoretical 
model  is  the  equation  of  state  (11-19)  relating  the  input  field  y  to 
the  output  field  x  as  a  function  of  the  inverse  of  the  effective  tem¬ 
perature  8s  • 

The  same  mechanism  is  well  known  to  be  responsible  for  the 
cooperative  atomic  interaction  which  results  in  the  dynamical  process 
of  superfluorescence  (23,31).  Our  model  also  connects  for  the  first 
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time  OB,  superfluorescence  of  the  superradiant  phase  transition  in 
thermodynamic  equilibrium  in  the  absence  of  applied  field  (32). 

A  completely  new  and  heretofore  unpredicted  phenomenon 
(29)  is  shown  in  Figure  4.  Here,  we  have  depicted  the  variation  of 
the  output  field  amplitude  x  as  a  function  of  the  inverse  tempera¬ 
ture  T  for  various  input  fields  yc  from  (III-5) .  These  curves  show 
bistability  and  reversible  hysteresis  in  the  output  field  x  with 
variation  of  the  inverse  temperature  r  for  fixed  input  field,  yc  for 
the  latter  greater  than  a  well-defined  threshold  value.  These  re¬ 
sults  could  have  important  applications  as  a  sensitive  FIR,  IR  or 
millimeter  wavelength  detector  with  frequency  upconversion.  That  is, 
an  FIR  or  IR  active  material  could  be  introduced  as  a  buffer  with 
the  optically  nonlinear  material  in  a  cavity.  The  FIR  or  IR  energy 
absorbed  by  the  buffer  can  be  kinetically  transferred  to  the  optical¬ 
ly  nonlinear  material  via  rotational-vibrational,  translational  ener¬ 
gy  transfer  to  alter  T  .  Switching  times  in  the  hysteresis  cycle 
could  be  on  the  order  of  the  thermal  fluctuation  times  for  the  system. 

We  also  predict  for  the  first  time  OB  in  a  small  volume 
without  mirrors.  For  the  case  without  mirrors,  the  mode  frequency 
0)  =  e  and  Y  =  c/2L  where  L  is  the  length  of  a  cylindrical  volume  of 
the  active  material.  The  conditions  for  a  cavity  are  that  the  mode 
frequency  w  is  the  cavity  frequency  and  need  not  equal  the  atomic 
frequency  e  ,  and  the  photon  escape  rate  Y  =  c/2Lt  where  T  is  the 
transmission  coefficient  for  the  cavity  mirrors.  Optical  bistability 
in  a  small  volume  without  mirrors  could  have  important  applications 
for  micron  size  OB  elements  imbedded  in  a  neutral  binder  to  consti¬ 
tute  a  two-dimensional  optical  digital  memory  bank  or  optical  digital 
imaging  surface.  This  certainly  would  alleviate  the  problem  of  hav¬ 
ing  to  produce  highly  uniform  reflective  surfaces  on  miniaturized  OB 
elements.  It  may  be  also  possible  to  constitute  a  three-dimensional 
optical  digital  memory  bank  of  relatively  compact  dimensions. 

Further  work  is  in  progress  (27,33)  to  calculate  the  rele¬ 
vant  switching  rates  in  terms  of  material  parameters  and  to  determine 
the  effects  of  atomic  line  shape  on  cooperative  OB.  We  are  also  in¬ 
vestigating  the  effects  of  multimode  coupling  to  the  Stark  shifted 
atomic  line  shape.  It  Is  anticipated  that  this  may  lead  to  the  es¬ 
tablishment  of  regenerative  oscillations  which  may  have  applications 
as  passive  high  intensity  laser  cw  to  pulse  conversion.  We  feel  that 
our  further  work  will  undoubtedly  lead  to  even  further  progress  and 
new  ideas. 
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